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ABSTRACT. We present a model of charmonium as two heavy quarks prop- 
agating classically in a weakly coupled quark-gluon plasma. The quarks 
interact via a static, color-dependent potential and also suffer collisions with 
the plasma particles. We calculate the radiation width of the color octet state 
(for fixed, classical qq separation) and find that it is long-lived provided a 
finite gluon mass is used to provide a threshold energy. 
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Theoretical studies of quarkonium in the quark-gluon plasma have gener- 
ally treated the qq bound state as a unit, a color-singlet meson propagating 
in the medium |TJ, |3|, [| , evolving slowly 0, |(| [7|, |8| , and colliding with the 
plasma's constituents || [L(| . Paradoxically, the interaction binding the qq is 



taken to be a static, screened Coulomb potential |TT|, 12 1, which implies that 



the density of plasma particles is high enough to allow construction of such a 
smooth potential — a weakly coupled plasma In this note we adopt 

the second point of view and construct a model accordingly. We consider [|15| 



a picture of quarkonium as two heavy particles propagating classically, sub- 
ject to a screened potential; these particles collide individually with plasma 
particles, changing momentum and color thereby. Such a model can be used 
to study in detail the diffusion of unbound cc pairs and the formation and 
destruction of charmonium in the plasma. 

For simplicity, we specify the potential in the plasma rest frame, and ne- 
glect the associated color-magnetic interaction. In writing the qq interaction 
we keep track of the color state c of the pair. The potential takes the form 

V c (r) = Cc— exp(-r/r D ) , (1) 
r 

where ( c = —4/3 and 1/6 for the singlet and octet states, respectively. We 
distinguish a priori three regimes in r (see Fig. 1): 

1. For large r the interaction is negligible and the system is color-blind — 
the two quarks diffuse independently in the plasma [16], and their color 



state (which is irrelevant to their motion) is determined statistically. 

2. For intermediate r one must keep track of the color state in detail, 
and propagate the quarks according to the potential appropriate to the 
instantaneous color state. 

3. For small r the large gap between singlet and octet might induce a large 
transition rate due to gluon radiation. If the width r 8 of the octet state 
is larger than the gap AV^, the octet ceases to exist as an excited state 
of the pair, vanishing into the (qq) S ingiet + 9 continuum. 

The boundaries of these regions must be determined by explicit calculation 
of transition rates. 

Before presenting the formal kinetic theory, we address the question of 
the existence of Region 3. Transitions between the color states are caused 
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by collisions of the individual heavy particles with plasma particles, and by 
radiation and absorption of gluons. Only radiation contributes to the octet 
state's width in the lowest order in a s . The radiating quark moves in the 
potential field of its partner, which is otherwise a spectator to the process 
(see Fig. 2). For an initial quark momentum p relative to the plasma, the 
radiation rate is given by 



X 



5(k + E p/ + V 1 (r)-E p - V 8 (r))j2\M\ 2 . (2) 



Here g(k) is the statistical factor for the outgoing gluon; the matrix element 
squared is averaged over initial spin and summed over final spins of the 
gluon and the participating quark. The summation over colors must be done 
carefully because of the condition that the initial state be an octet and the 
final state, a singlet. We define a basis of color state vectors \ij), where 
i, j = 1, 2, 3 are the color indices of the quark and antiquark, respectively. 
Then the singlet state is 

|1> = ^EI«> ( 3 ) 
while a member of the octet is[] 

|8r) = >/2X;W> • (4) 

u 

The Feynman matrix element then contains the color factor 

V6 

where the color operator T a acts on either the quark or the antiquark state 
only. Squaring the matrix element, summing over the color a of the emitted 
gluon, and averaging over r gives 



3 We normalize the color matrices so that Tr T a T b = \5 ab - 
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The result is 

J2\M\ 2 = 8 -^(2m^p.p>) . (7) 

Of course, either the quark or the antiquark may radiate; if their momenta 
are p\ and P2, then the total radiation rate from the pair is r 8 = r(pi)+r(p 2 ). 

Calculating T{p) from (||) gives the nonsensical result that T{p) 7^ in 
the limit AV — > 0. In this limit 

The source of the difficulty is that the Bose-Einstein factor contributes 1/k to 



counteract the vanishing phase space JL7|] . It is evidently a bad approximation 



to take the radiated gluon to be massless; the mass induced by the plasma 
must be taken into account. The simplest thing to do is to take (0) and to 



replace k in the measure, in g(k), and in the energy S function by J k 2 + m 2 . 
This makes Y vanish for gap values below the threshold AV shown 
in Fig. 3. A more correct calculation would replace the energy 5 function 
by the one-loop spectral density of transverse gluons (18], [L9] , and include as 



well the radiation of longitudinal plasmons, similarly treated. 

The qualitative result indicated by Fig. 3 is that there is no distinct small- 
r regime as we supposed. The octet state is meaningful even when AV is 
large. Of course, T 8 will be increased by the effects of collisions, and so it is 
still possible that the octet disappears at some value of r. 

The radiation and collision rates are ingredients in the Boltzmann-Vlasov 
equation which governs the evolution of the qq distribution. This distribution 
is defined as a function / c (xi, x 2 , p 1; p 2 ; t) on the two-particle phase space, 
carrying as well the index c denoting the color state of the pair. It satisfies 
the transport equation 

JMf) • (9) 

\ / collision 

The left-hand side is the convective derivative in phase space, 



Dt J 



d pi d p 2 d „ . , / d d 

at m ax! m ax 2 Vc'Pi 
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where F c (r) = — VV c (r) is the color force. The right-hand side of ([|) rep- 
resents gain and loss through collisions and radiation. In the small- and 
intermediate-distance regimes, we distinguish collisions which change the 
color state from those which do not. The former (like the radiation pro- 
cess) carry momentum-transfer of at least AV(r), and hence are hard colli- 
sions; the latter may carry very small momentum transfer, and hence may 
be approximated as soft. Thus we write 



'dfc\ 



— -Rsoft + -Rhard ■ (H) 

collision 



We can approximate the soft collision term by Fokker-Planck terms [16 



= - £ ' P-WJ + aS^W /J + (1 - 2) • (12) 

G c and are drag and diffusion coefficients, calculated in the Landau ap- 
proximation as moments of transition rates [^, [16[]. Here, as for the radiation 
rate calculated above, one has to be careful in projecting out the proper color 
channel. This makes G c and dependent on c, and implies that only a 
subset of the scattering diagrams contributes to a given color channel. 

We have calculated the drag and diffusion coefficients in the two color 
channels [pl|. These can come only from soft collision processes, since ra- 
diation necessarily flips color. The octet coefficients are similar to those for 



the color-blind scattering which a single quark undergoes [|6[], reduced by 
roughly a factor of 2-3 because of the restriction of the final states. The 
coefficients for the singlet state are smaller by an additional factor of 100, 
again because of the reduced number of final states, and also because of the 
smaller set of diagrams which contribute. 

The hard collision terms cannot be approximated, and must be written 
as full-blown collision integrals. We have the loss and gain terms, 

-Rhard = - [r c (pi) + r c (p 2 )] / c (x , p a ; t) 

+ J d 3 k w(px - k -> px ; c -> c) / E (x a , p x - k, p 2 ; t) 

+ J d 3 kw(p 2 - k -> p 2 ; c -> c) /c(x a ,pi,p 2 - k;i) . (13) 

Here T c (p) is the total rate of radiation- and collision-induced transitions 
for each quark, along the lines of T defined in (0). The functions w are 
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differential transition rates into the phase space element at (x a , p a ) from 
other elements; since these are hard collisions, they involve perforce a color 
flip into c from the complementary state c, and depend on position through 
the value of AV(r). 

Our result for the radiative width of the octet state gives a lifetime of 
several fm/c for reasonable qq separations. If the collision rates do not affect 
this result strongly, then the evolution of the pair in phase space requires 
detailed numerical solution of the Boltzmann equation. If, on the other 
hand, the collisional width is large, then it may be reasonable to assume that 
the singlet and octet states are in equilibrium. Then a simplified Boltzmann 
equation may be studied, in which f c oc exp(— V c /T). If the collisional width 
turns out to exceed the gap AV at some values of r, then the octet state 
ceases to exist and we are back to the consideration of three distinct regimes 
in r as described above. 
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Figure 1: The qq potential in color-singlet and color-octet channels, not 
drawn to scale. 




Figure 2: Feynman diagram for gluon emission from a quark, z, j, and a are 
color indices. The lower line represents the spectator antiquark. 



6 




Figure 3: Radiative width T 8 = T (pi) -\-T (p 2 ) of the octet state at pi = p 2 = 
as a function of the gap AV . We set a s = 0.6 and the gluon mass is 300 MeV. 



7 



